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Abstract 



Phase separation in mixtures of polymers and alkanes is investigated with the aid of 
a recently developed lattice model, based on the Born- Green- Y von (BGY) integral 
equation approach to fluids. The system-dependent parameters of the BGY lattice 
model for binary mixtures are deduced from those of its pure components, which in 
turn are determined from a comparison with experimental data. The lower critical 
solution temperatures (LCST's) for polyethylene in various n-alkanes were predicted 
from the BGY lattice model and compared with experimental data. While the model 
underestimates the LCST's for the smaller alkanes it reproduces the experimental 
values very well for the higher alkanes (decane through tridecane) . The effect of the 
chain length of the polymer on the solubility is investigated for the case of decane as 
a solvent. 

Introduction 

In many polymer /solvent systems, phase separation occurs at a lower critical solution 
temperature (LCST) 0]. For temperatures smaller than the LCST the polymer is sol- 
uble while a miscibility gap opens above the LCST which is closed again at an upper 
critical solution temperature (UCST). These lower critical solution temperatures are 
difficult to predict accurately from theory 0, Q . In this work, we employ a recently 
developed lattice model HI, H, ^ H based on the Born-Green- Yvon (BGY) integral 
equation approach to fluids to investigate the miscibility of polyethylene in alkanes. 
The BGY lattice model provides a compact expression for the configurational part of 
the Helmholtz free-energy density which allows the calculation of the relevant ther- 
modynamic properties of a mixture. In addition to the system-dependent parameters 
of the pure components, the lattice model for a binary mixture contains an energy 
parameter describing interactions between unlike molecules. In this work the inter- 
action energy is approximated by the geometric mean of the energy parameters of 
the pure components, allowing the thermodynamic properties of the mixture to be 
predicted from those of the pure components. In the following sections we describe 
how the BGY lattice model is used to predict the LCST's and corresponding coexis- 
tence curves for polyethylene in various alkanes, and then compare our results with 
experimental data. 
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Born-Green- Yvon lattice model for fluids 



In the BGY lattice model, a molecule of species i is assumed to occupy contiguous 
sites on a lattice having coordination number z and a total of Nq lattice sites. As 
in earlier work 0, |], we set z = 6. The site-fraction 0j is defined as 0j = ViNi/No, 
where Ni is the number of molecules of component i while 0h = Ny^/Nq is the fraction 
of empty sites or holes, where A^h is the number of holes. Each molecule of species 
i has qiZ = ri{z — 2) + 2 interaction sites, leading to the definition of concentration 
variables = qiNi/{N-^ + X]j Ij^j) and = A'h/(iVh + Hj Qj^j) which account for 
the nearest-neighbor connectivity of the molecules. The interaction energy associated 
with nonbonded nearest neighbors of the same species is en, while eij corresponds to 
interactions between unlike nearest neighbor segments. Using the Born- Green- Yvon 
integral equation approach, Lipson [^, § derived an expression for the dimensionless 
configurational Helmholtz free energy per lattice site a = (3A/Nq of an M-component 
mixture 
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where /3 = 1/(A;bT), T is the temperature and /cb is Boltzmann's constant. Denoting 
the total volume of the lattice by = vNq, where v is the volume per site, we obtain 
the Helmholtz free-energy density A/VT 



A/VT = k-Qa/v + Oq, 



(2) 



where is a caloric background which is not discussed further in this work. All 
thermodynamic properties of a mixture can now be derived from the thermodynamic 
relation 
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where U is the internal energy of the system, and where pj = Nj/V = (pj/vjV and 
fij are the number density and chemical potential of component j, respectively. The 
pressure P and the chemical potentials pi, for example, are given by 
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It now becomes convenient to focus on molar quantities, so that pj becomes the 
molar density of component j, v denotes the volume of a mole of lattice sites, the 
energy parameters e are measured in J/mol, and k-Q is replaced by the universal gas 
constant R. 





Application to one-component fluids 

For the case of a one-component fluid, equations are evaluated with M = 1, 

so that we can drop the subscripts other than h. There are three system-dependent 
parameters which we determine from a comparison with experimental data: the num- 
ber r of lattice sites occupied by a single molecule, the volume v per mole of sites, 
and the interaction energy e. 

To obtain parameters for the ?7,-alkanes from pentane through tridecane we con- 
sidered coexistence densities. The critical temperature and the critical density Pc 
of a one-component fluid are determined by the following two conditions: 

= 0, (5) 

For temperatures T smaller than Tc the coexistence densities pvap and puq can be 
found from the conditions that pressure and chemical potential are equal in the vapor 
and the liquid phase, i.e. 

P{T, Pvap) = P{T, pliq), p(T, Pvap) = p(T, pUq) . (6) 

We determined the parameters r, v, and e for a number of ra-alkanes between methane 
and octane by comparing calculated liquid and vapor densities for temperatures be- 
tween 0.5Tc and 0.9Tc (to avoid the critical region) with values tabulated by Var- 
gaftik [0 . We then correlated the parameters with the chain length n of the alkanes, 
where we typically excluded the results for the lowest alkanes. The resulting expres- 
sions for the parameters as a function of chain length are given by 

rt;/(L/mol) = 0.016n + 0.015, (7) 
re/(J/mol) = -2150n - 3150, (8) 

^ = -0.877539 + 1.512447Vn + 1.4651 ( ^ + 0.58259 ) . (9) 
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Figure 1: The coexisting vapor and liquid densities for n-hexane. The symbols indi- 
cate values from Ref. 0, the solid line values calculated from Eqs. @ with system- 
dependent parameters given in Table |1]. 

The first of these equations, Eq. (|^, reflects the fact that the hard core volume of the 
alkanes is expected to grow linearly with chain length, while Eq. (P) conveys a linear 
relationship between the interaction energy per mole and the chain length n. Due to 
end effects, neither relation is expected to hold for very small n. The last equation, 
Eq. @, is motivated by the variation of the experimental |]10[ and calculated critical 



temperatures of the alkanes with chain length. In Table |1| we present the values of the 
system-dependent parameters for pentane through tridecane according to Eqs. (|^)- 
(H). Calculated values for the coexistence densities of hexane along with values from 
Ref. are presented in Fig. while percent deviations between calculated and 
tabulated values for all alkanes considered in this work are presented in Fig. |^. The 
highest deviations seen in Fig. |^ are for vapor densities and are, in part, a result of 
the correlation of the parameters with chain length. The two flgures show that our 
model gives a reasonable description of the coexistence densities outside the critical 
region. 

System-dependent parameters for polyethylene (PE) were determined from a com- 
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Figure 2: Percent deviations between values for the coexistence densities as provided 
in Ref. |^ and as calculated from Eqs. with system-dependent parameters given 
in Table |l|. 



Table 1: System-dependent parameters for the pure components. 



name 


molmass 


r 


V (mL/mol) 


e (J/mol) 


LCST(K) 


n-Pentane 


72.146 


9.505229 


9.994499 


-1462.353 


260.4 


n-Hexane 


86.172 


10.35107 


10.72353 


-1550.565 


371.2 


n-Heptane 


100.198 


11.22030 


11.31877 


-1622.060 


435.9 


n- Octane 


114.224 


12.10191 


11.81632 


-1681.552 


484.0 


n-Nonane 


128.250 


12.99011 


12.24009 


-1732.088 


522.8 


n-Decane 


142.276 


13.88162 


12.60660 


-1775.729 


555.5 


n-Undecane 


156.302 


14.77453 


12.92766 


-1813.933 


583.7 


n-Dodecane 


170.328 


15.66765 


13.21193 


-1847.756 


608.6 


n-Tridecane 


184.354 


16.56028 


13.46595 


-1877.987 


630.8 


PE 


140262.0 


15890.72 


9.621464 


-1988.320 




PE 


85000.0 


9629.90 


9.621934 


-1988.458 




PE 


50000.0 


5663.81 


9.623040 


-1988.729 
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Figure 3: Percent deviations between experimental [|TT], |l2l and calculated (cf. 



Eq. (^)) values of the densities of polyethylene at given temperature and pressure. 



parison with experimental PVT data [TT|, [12| in the liquid region for temperatures 



between 420K and 560K and pressures between lOMPa and 200MPa. When the 
density of a long chain polymer is given as a mass density, rather than a number 



density, the pressure is expected to be independent of the molecular mass |12[ of 
the polymer. We determined system-dependent parameters for PE from a compar- 
ison between experimental and calculated pressures using three different molecular 
masses for the polymer. The representation of the experimental data was virtually 
identical in the three cases, and the resulting values for the parameters, as presented 
in Table |l|, show that for long chain polymers e and v are indeed independent of the 
molecular mass, while r is proportional to it. A comparison between calculated and 
experimental values for the density at given temperature and pressure is presented 
in Fig ^. The deviations are generally smaller than 0.3% but reach up to 1% for the 
lowest pressures. 



Phase separation in binary solutions 



To discuss phase separation in binary mixtures, it is convenient to introduce the total 
molar density p = pi + P2 and the mole fraction x = P2/ p of component 2. The 
thermodynamic field conjugate to the mole fraction x is the difference of the chemical 
potentials p = p2 — Pi- In terms of these variables, the conditions for a point Tc, Pc, 
Xc on the critical line of a binary mixture take the form []T3| : 





T,P \""^/r,p \^t^/T,xl / T,x 



0. (11) 



T,P 



The first of the two equations defines the spinodal while the second one implies that, 
for given pressure, the critical point is an extremum of the spinodal in the T-x plane. 
Hence, lower and upper critical solution temperatures (LCST's and UCST's) corre- 
spond to the minimum and maximum of their spinodals, respectively. For tempera- 
tures LCST < T < UCST, the mixture phase separates into two phases of different 
mole fractions xj and xjj. For a given pressure P, the coexisting phases satisfy 

p{T, P, xi) = p{T, P, xii), pi{T, P, xi) = /ii(T, P, xn). (12) 

For typical solutions of polymers in alkanes, the mole fractions x are very small 
numbers so that it is more convenient to employ mass fractions c: 

c = M2x/(M2X + Mi(l-x)), (13) 

where Mi and M2 are the molecular masses of component 1 and 2, respectively. 

When the BGY lattice model is applied to binary mixtures, Eqns. (|l|)-(^) are 
evaluated with M = 2. The total molar density p = pi + P2 and the mole fraction 
X = P2/P of the mixture are related to the site fractions 0i and 02 by 

, = if*+!^), i = l + (14) 

V \ri r2 J X ri 02 

In addition to the energy parameters en = ei and 622 = ^2, and the numbers of 
segments ri and r2, which are determined directly by the pure components of the 
mixture, values need to be supplied for the interaction energy eu and the volume v of 
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a mole of lattice sites. For the energy parameter eu we employ Berthelot's geometrical 
mean combining rule [13]: 



ei2 = -y/^- (15) 

To find a value for v, we take advantage of the fact that the BGY model for the pure 
polymer depends more strongly on the product rv than on the values of the number 
of segments r and the volume v separately. Hence, for a mixture of an alkane and 
polyethylene, we assign the value Vi of the alkane to the mixture, and rescale the 
segment length of the polymer, i.e. 

v = vi, {r2V2)/vi. (16) 

When parameters v and r2, given by Eq. (|16D, are used to compare experimental and 
calculated PVT data for pure polyethylene over the whole temperature and pressure 
range, it is seen that the new parameters diminish the overall agreement between 
experiment and theory. However, for the lowest pressures the description of the PVT 
surface is improved over that using the original parameters, and becomes excellent 
for the highest alkanes. This is important since we are going to investigate phase 
separation of polyethylene solutions at atmospheric pressure. 

To determine the lower critical solution temperatures for mixtures of polyethylene 
and the various alkanes discussed in the previous section, we employed the parameters 
presented in Table |l|, where the highest molecular-mass parameters were used for 
polyethylene, together with Eqs. ([15[) and (piBI) . We then solved Eqn. ( |T0|) numerically 



to find the spinodal curve for atmospheric pressure P = O.lMPa. The LOST was 
identified as the minimum of the spinodal curve as a function of x. All calculated 
LCST's are included in Table |I], and the values for the n-alkanes between pentane and 
tridecane are compared with experimental values [0, ^, |T^, in Fig. ^ These results 



show that while the BGY lattice model with the parameters of Table |l] somewhat 
underestimates the LCST's for the shorter chains, it gives a very good description of 
the data for the higher alkanes. 

The influence of the chain length of the polymer on its solubility is illustrated in 
Fig |[ The symbols indicate experimental data by Kodama and Swinton 0] for three 
polyethylene samples of different molecular masses in decane. We calculated the coex- 
istence curves for polyethylene in decane for three comparable molecular masses (see 
Table |I|) using the coexistence conditions Eqs. ([T^p. Shown in Fig. |^, as a function of 
the mass fraction c, are the differences between the coexistence temperatures and the 
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Figure 4: Lower critical solution temperatures for polyethylene in n-alkanes at P = 
O.lMPa. The symbols represent experimental data for Marlex-50 by Orwoll and 
Flory [0 (open circles), for Kodama and Swinton's longest chain polymer PE 
Type 2 (open squares), extrapolated values to M — > oo by Hamada et al. (open 
diamonds), and values collected by Charlet and Delmas [|1^] (open triangles). The 
filled diamonds on the solid line indicate values calculated with the aid of Eqns. (|l^ 
and ( p!Ol) and the values presented in Table |l] together with Eqns. (|T5|) and (|16|) . 
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Figure 5: Coexistence curves for polyethylene in n-decane at P = O.lMPa as a func- 
tion of the mass fraction c. The symbols correspond to experimental data by Kodama 
and Swinton the lines represent values calculated according to Eqn. (0) with the 
aid of the parameters in Table |I| and Eqns. ( ]TB| ) and (|T^). ^ For the experimental 
data, the LCST (T = 557.6K) of PE Type 2 has been subtracted from the coexistence 
temperatures, while the predicted LCST (T = 555. 7K) of PE with molecular mass 
140262 has been subtracted from the calculated values. 
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calculated and experimental values of the LCST of the longest chain, respectively. 
Since the experimental data in the paper by Kodama and Swinton 0] were not tabu- 
lated but presented in a graph the comparison is not exact. It appears, however, that 
the effect of changing molecular mass on the coexistence temperatures of polyethylene 
in decane is predicted well by the BGY lattice model. 

Conclusion 

With the aid of the system-dependent parameters for the pure components as well 
as the geometric mean approximation for the interaction energy, the BGY lattice 
model has been shown to predict LCST's for solutions of polyethylene in alkanes in 
very good agreement with experimental data. The theory also captures the effect of 
changing polymer molecular mass on solubility. In current work, we are investigating 
the influence on miscibility of polymer and solvent architecture. 
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